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Groundstate properties of the BCS-Bose Einstein crossover in a d x 2_ y 2 wave 

superconductor 
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We examine the groundstate properties of the crossover between BCS superconductivity and Bose 
Einstein condensation within a model which exhibits d, x i_ y ? pairing symmetry. We compare results 
for zero temperature with known features of s-wave systems and show that bosonic degrees of 
freedom are likely to emerge only in the dilute limit. We relate the change in properties of the 
groundstate as a function of density to the effect of the exclusion principle and show that the one 
particle distribution function undergoes a significant change when bosonic behavior appears. 
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I. INTRODUCTION 

The existence of a crossover from BCS supercon- 
ductivity to Bose Einstein (BE) condensation of pre- 
formed pairs has received increasing attention in the 
literature .over the last few years. Originally discussed 
by Eaglesu within the context of pairing in thin film 
semi-conductors, -and considerably expanded upon by the 
works of LeggettEl and Nozieres and Schmitt-RinkH, the 
latest resurgence in interest has been due to its possible 
application to the understanding of the phase diagram of 
high temperature cuprate superconductors. 

In particular, the appearance of a pseudo-gap ipj-both 
charge and spin excitations of the normal stateQcl has 
led to a number of suggestions that pairing correlations 
well above T c may occur in these systems. Whether 
these-,c«rrelations are in the bosonic form of pre-formcd 
pairsatJJ, or pair resonances originating from the interme- 
diate regime of a BCS-BEresossoveiOO, or simply from 
classical phase fluctuations^ is still however a controver- 
sial issue. 

Furthermore, now that there is a large body of evidence 
to suggest that pairmgria the cuprates is predominantly 
d x 2_ y 2 in char act enljiiZl, there is a clear need to un- 
derstand the properties of the BCS-BE crossover within 
the context of this type of pairing symmetry. Although 
there have been some attempts to discuss the effect of 
d x 2_ y 2 pairing on pseudo-gap formation above T c in the 
cupratestS, there has been little discussion on the sys- 
tematic groundstate properties of the BCS-BE crossover 
in the d x 2_ y 2 channel. This is an important issue, since 
pairs with this symmetry cannot contract to point like 
bosons and so accordingly one expects that there will 
be severe consequences for the properties of the ground- 
state of this type of system. This is of direct importance 
in terms of the validity of the crossover scenario for the 
cuprates, whilst also being of general interest in under- 
standing macroscopic pairing in higher angular momen- 
tum channels. 

In this paper, we consider the BCS-BE crossover at 



zero temperature as a function of both coupling strength 
and carrier density, within a two-dimensional toy model 
which has a d x 2_ y 2 pairing instability. The study at zero 
temperature is well controlled by use of the BCS_varia- 
tional wavefunction (which contains the BE limitcl) and 
allows us to establish at the two body level the qualitative 
groundstate properties of the system upon which future 
studies may be based. In particular we show that the 
groundstate properties of the system are severely modi- 
fied from the s-wave picture, where a smooth crossover 
exists for all densitiesBu. For the d x 2_ y 2 case the ef- 
fect of the exclusion principle as the density of carriers 
is increased results in a suppression of the emergence of 
bosonic degrees of freedom for moderately large densities 
and strong coupling, where the system instead remains 
fermionic (BCS like). We find that only in the dilute 
limit is a crossover possible, and that when this occurs 
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FIG. 1. The crossover between fermionic and bosonic de- 
grees of freedom as a function of carrier density and coupling 
strength (as defined in text). The left boundary is deter- 
mined by the onset of a finite gap amplitude, and the right 
boundary by when the chemical potential falls below the band 
minimum. The solid lines are a guide for the eye only. 
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the single particle distribution function undergoes a rad- 
ical change. We summarize most of our results in Fig. [I], 
which shows when bosonic degrees of freedom can emerge 
as a function of both density and coupling strength. 



II. TOY MODEL 



We introduce as our toy model a 'reduced' Hamilto- 
nian in the BCS sense which describes an effective two- 
particledxtteraction in real space and in the singlet pairing 
channeO. Specifically it is 



H = ^ ~t( C lcr C jcr + h.C.) — /X* rii a 
<ij>(T ia 

i <ij> 

where first and second terms describe nearest neighbor 
hopping on a two-dimensional square lattice with chemi- 
cal potential fj,* , and the third and fourth terms describe 
an effective two-body pairing interaction in real space. In 
particular, the third term represents the repulsive part of 
the effective interaction while the last term provides, an 
attractive interaction for nearest neighbor particlesel 

By expressing the superconducting gap in terms of its 
various symmetry components and introducing the BCS 
variational wavefunction \^) = Jlk( M k + ^kcj c |cl_ k |)|0}, 
the zero temperature gap equation in the d x i_ y i channel 
is given by, 



V (cos k x — cos ky) 2 



1 = J-V 

2M k K£(k) - M) 2 + A 2 (cosfc a: - cosfcj,) 2 



1/2 



, (2) 



where £(k) = —2tr]k is the nearest neighbor tight- 
binding energy dispersion, the geometric factor ?/k = 
cos k x + cos k y and the effective chemical potential [i = 
/i* — n(W/2 — 2V). The amplitude of the d-wave gap is 
denoted by A. n 

As Eagles first pointed outfil, any deviation from weak 
coupling requires a self consistent solution of both gap 
and number equations, since the BCS approximation of 
the chemical potential being equal to its value in the nor- 
mal state can no longer be reasonably justified. Specifi- 
cally, the number equation for the chemical potential /i, 
which defines the particle density n = N/M is , 



n - 1 = 



-E 



-(C(k)-M) 



[(^(k)- A1 )2 + A2(cOSfc :c 



cos ky) 2 ] 1 / 2 
(3) 



We have solved Eqns. (g) and (J3|) self consistently at 
a given density by numerical integration. In this theory, 
pairing takes place over the entire Brillouin zone (BZ). 
The zone edge acts as a natural boundary or momentum 
cut-off which avoids the renormalization methods used to 
remove ultraviolet divergences in continuum model tr,&at- 
ments of strong coupling within the BCS frameworkEil. 



III. CROSSOVER AND ANALYSIS 

The inset of Fig. || shows the gap parameter A for 
densities ranging from the dilute limit to the almost half- 
filled case. Taking the onset of a finite gap parameter as 
the signal for the manifestation of the superconducting 
state, an increase in the particle density clearly favors 
the emergence of a d x i_ y i paired groundstate. On the 
other hand, as the system becomes more dilute, there is 
a need for stronger coupling to induce pairing. 
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FIG. 2. The dependence of /i on coupling strength for var- 
ious densities. Inset: The gap A as a function of coupling 
strength (dashed lines and data points). The solid lines are 
the asymptotic behavior derived in the text. 

In the dilute limit (of relevance to the underdoped 
cuprates, where the density of carriers is proportional to 
the doped hole concentration x), the occurrence of d x 2_ y 2 
bound states for the two-particle problem on an empty 
lattice is for our system, equivalent—to that problem for 
the t — J model. Kagan and RiceE3 have shown that a 
d x 2_ y 2 bound state will not occur unless the coupling J 
(equivalent to V in out. case) is greater than V c /4t ~ 1.8. 
Also, Randeria et alB2 have shown that in two dimen- 
sions, a necessary and sufficient condition for a dilute 
many-body s-wave Cooper instability to occur is that an 
s-wave bound state exists for the corresponding two-body 
problem on the empty lattice. Importantly, in the con- 
text of d x 2_ y 2 symmetry, they have also shown that such 
a condition does not exist for higher angular momentum 
pairing. 

In our study, we find that the onset of pairing occurs 
for progressively weaker coupling as the density is in- 
creased, contrary to an s-wave system. In the d x 2_ y 2 
channel at a coupling strength less than V c /4t, the dilute 
system gains more kinetic energy compared to pairing 
energy and thus has a total condensation energy which 
is positive, whereas for the more dense system the ki- 
netic energy of the carriers is on average less effected 
and a pairing instability occurs. This steady evolution 
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from the extreme dilute limit where pairing does not oc- 
cur until V/At > 1.8, to the near half filled case where 
superconductivity can manifest itself at a much weaker 
coupling than that required for a two body bound state, 
is indicated by the boundary on the left in Fig El The 
n = point corresponds to the critical coupling strength 
V c discussed above and calculated initially by Kagan and 
RiceM 

In Fig. we show [i as found in conjunction with 
the solutions for the gap A at various densities and as a 
function of the coupling strength. The horizontal line at 
fi/At = —1.0 represents the bottom of the tight binding 
band. For weak coupling, it is well known that /i in the 
superconducting phase is given roughly by the Fermi en- 
ergy of the normal state and this can be seen in the figure. 
However, at large doping, jj, shows little deviation from 
its normal state value over a large variation in the the 
coupling strength. This is to be contrasted with the low 
density results, which show a relatively rapid deviation 
from weak coupling behavior as V is increased. 

In general, bosonic degrees of freedom can be expected 
to emerge once the chemical potential of the many-body 
groundstate slips below the band minimum in a tight 
binding system, or below zero in a continuum paodel. 
For an s-wave system, Nozieres and Schmitt-RinkfcJ were 
able to show that a crossover from fermionic supercon- 
ductivity to bosonic degrees o freedom can occur for all 
densities as the coupling strength is increased. For the d- 
wave system considered here, this is not the case. Bosonic 
degrees of freedom can only emerge in the dilute regime, 
while for large densities, the system behaves more like a 
weak coupling superconductor with a value of the chem- 
ical potential comparable to that of the normal state. 
These results are expressed by the boundary on the right 
in Fig. [lj Given that the BCS wavefunction only takes 
two-body correlations into account, we expect that the 
suppression of bosonic degrees of freedom would be in- 
creased in a more sophisticated treatment which would 
include a repulsive interaction between fermion pairsEj. 

Further insight into this feature of a d-wave system 
can be gained by examining the limit of infinite coupling 
strength. In the case V — ► oo, the kinetic term (and thus 
any Fermi surface geometry) becomes negligible and the 
asymptotic behavior of the gap and chemical potential 
can be shown from Eqns. (|J) and (||) to have the form 
A — > r yV/2 and — > 7V r (n — l)/2a, implying that in 
the infinite coupling limit, /x/A — > (n — l)/a. Here the 
parameter a is given by the solution to 



(4) 



while 7 is defined by 



cos ky) 2 



({n-l]/a 



,211/2' 



7 = 



M ^ 



(cosk x — cos ky )" 2 



[(cosk x — cos k y ) 2 + ([n — l]/a) : 



a/2- 



(5) 



We have indicated the asymptotic behavior of A at 
various densities in the inset of Fig. H by the solid 



lines. As one increases the density, the convergence to 
the asymptotic behavior is poor for large densities, indi- 
cating that the fi is still of the order of the band energy. 
In the strong coupling regime, if true bosonic character- 
istics emerge then one expects that /i should simply re- 
duce to the binding energy per particle for the 'diatomic' 
electron molecule. It can readily be shown that only in 
the dilute limit does the ratio 7/a approach unity. This 
indeed leads to the result = —V/2, exactly one half 
of the binding energy-jfor the two body problem in the 
strong coupling limited. On the other hand, at half fill- 
ing (n = 1) fi clearly remains zero even in the limit of 
infinite coupling. Remarkably, the groundstate remains a 
fermionic superconductor for all coupling strengths and 
has characteristics equivalent to a weak coupling BCS 
system. For densities in between these two limits, an 
increase in coupling will eventually reduce the chemical 
potential below the bottom of the band, however it will 
always be somewhat greater than the binding energy per 
particle of the two-body case. 

This general dependence on the density for the BCS- 
BE crossover in ci-wave systems can be viewed as a man- 
ifestation of the effect of the exclusion principle. Due to 
the symmetry of the pairs, they cannot contract in real 
space to point bosons, but must always retain a finite 
spatial extent. As the density increases, the overlap of 
the pair wavefunctions exerts its influence through the 
exclusion principle contributing a positive energy to the 
system. At half filling, this overlap prevents the system 
from crossing over to a system displaying bosonic quali- 
ties even in the infinite coupling limit. 

This interpretation can be further supported by a cal- 
culation of the average pair coherence length. We can 
define this length £0 through the expectation value of the 
quantity $ = (F k | - V£ |F k )/(F k |F k ) A«here F k = u*t; k 
plays the role of the pair wavefunctionE3. 
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FIG. 3. The coherence length in the strong coupling limit 
as a function of density. The lattice spacing has been set to 
unity. 
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In Fig. |3j we show the behavior of £o in the infinite 
coupling limit as a function of density. It is clear that a 
shrinking of the pair size in real space to a compactified 
boson spread out over nearest neighbor sites can only 
occur in the dilute limit, where the exclusion principle 
becomes irrelevant. For higher density, the overlap of 
the pair wavefunctions has the effect of increasing the 
correlation between particles, which then increases the 
average pair size in order that the total energy of the 
system can be minimized. 

Lastly, we examine the single particle distribution 
function function of the coupling strength. We 

find that in weak coupling, the system behaves as a typ- 
ical BCS system, with riko- resembling a step function 
up to a point close to the normal state Fermi surface, 
where upon it becomes smeared over a small region about 
fx = Ep. The degree of smearing increases with the cou- 
pling until the chemical potential drops below the bot- 
tom of the band. At this point, the behavior of the single 
particle distribution function radically changes. 

In Fig. |J, we compare contour plots of for n = 0.3 
throughout the BZ for the weak coupling case V/At = 1, 
where \x is well approximated by its normal state value, 
and the stronger coupling case V/At = 4 where the chem- 
ical potential falls just below the bottom of the band. 
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FIG. 4. Contour plots within the first BZ of for 
V/At = 1 (top) and V/U = 4 (bottom) and n = 0.3. The 
brighter the region the larger the value of nk CT . 

In the weak coupling case, the Fermi surface of a tight 
binding band with a small density per unit volume can be 



clearly seen by the bright region in the middle of the plot. 
On the other hand, when fj, drops below the bottom of the 
band and bosonic degrees of freedom emerge, nko- under- 
goes a redistribution within the BZ. For V/4t = 4, the 
probability for occupation of highest momentum states 
is now found in the regions around (±7r, 0) and (0, ±7r). 
For d x 2_ y 2 pairing, this change in behavior of the single 
particle distribution function as the chemical poten- 
tial falls below the minimum of the tight binding band is 
an interesting feature. For s-wave systems, nko- becomes 
a constant for strong coupling, representing the Fourier 
transform of a point internal wavefunction. Accordingly 
we can interpret the new structure for n\^ a as the d-wave 
version of a local pair. 

One may speculate what the effect of this new struc- 
ture for riko- may be at finite temperatures. Above T c , if 
a pseudo-gap in the normal state excitation spectrum can 
arise from pairing fluctuations in the crossover regimelia, 
then one would expect these fluctuations to occur pre- 
dominantly in the region of the BZ which has the largest 
probability of occupation by pairs. On the one hand, it 
is interesting to note from the bottom plot in Fig. ^ that 
these regions correspond to the angular dependence of 
the pseudo-gap in the underdoped cupratesErcl. However, 
it is still unknown whether these materials correspond to 
such a regime. 

IV. SUMMARY 

In summary, we find that for d x 2_ y 2 pairing, only in the 
dilute limit is it likely that a BCS-BE crossover can occur, 
while it is possible at any density for s-wave systems. 
If bosonic behavior does emerge, the d x 2_ y 2 symmetry 
causes the single particle distribution function to undergo 
a radical redistribution. 
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